Optimal Stochastic Enhancement of Photoionization 
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The effect of noise on the nonlinear photoionization of an atom due to a femtosecond pulse is 
investigated in the framework of the stochastic Schrodinger equation. A modest amount of white 
noise results in an enhancement of the net ionization yield by several orders of magnitude, giving 
rise to a form of quantum stochastic resonance. We demonstrate that this effect is preserved if the 
white noise is replaced by broadband chaotic light. 
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The interplay between noise and the nonlinear re- 
sponse of a physical system described by classical me- 
chanics has led to intriguing effects. A paradigm is 
stochastic resonance (SR), whereby a nonlinear system 
responds more strongly to coherent subthreshold driv- 
ing in the presence of an appropriate amount of noise 
[E S 9 • Despite prominent demonstrations of classical 
SR in ample variety, only a few quantum mechanical ex- 
amples have been studied possibly due to subtle features 
of the quantum evolution [3| . The concept of SR in the 
quantum domain was originally suggested in the context 
of two-level conductance switching of driven mesoscopic 
metals 4] . This proposition sparked a number of theoret- 
ical studies on quantum stochastic resonance in double- 
quantum well structures [5| , the so-called incoherent and 
coherent SR in driven dissipative spin-boson systems 6 1 , 
in the bistable photon-field dynamics in micromaser 7], 
and in the electron shelving in a single ion Q. Experi- 
ments using NMR spectroscopy of water and a very 
recent one in the radiofrequency controlled nanomechan- 
ical beam oscillator |10|] , have established the properties 
of quantum SR in two- level systems (TLS). These stud- 
ies are mostly restricted to the quantum analog of the 
classical double- well dynamics. Yet, they have provided 
valuable insight into the noise-induced dynamics of quan- 
tum systems. 

Following quite a different path of research, there is 
a rapidly growing activity in the general area of control- 
ling quantum phenomena llj. A common approach to 
exercise the control exploits the non-perturbative inter- 
play between a purposefully designed optical field from 
a laser and the target quantum system, such as an atom 
or molecule [12J. A fundamental control goal is the ma- 
nipulation of atomic and molecular interaction to steer 
the quantum system towards a desired state [ill [T3 | . An 
accurate knowledge of the effect of noise on quantum sys- 
tems would be very helpful to achieve full control. Their 
response to noise has been rarely studied so far [l3] . One 
may even wonder whether the presence of noise offers 
new possibilities of quantum control. 

Here, we will demonstrate the existence of a stochas- 
tic resonance-like effect in a generic quantum situation 



beyond the two-level systems. For this purpose, we con- 
sider a quantum system having a finite binding potential 
with mixed, discrete and continuous spectrum which is 
coupled to two external forces. First, a nonresonant co- 
herent optical driving and, second an incoherent pertur- 
bation which may result from some form of environment. 
These situations are sufficiently general to be achieved in 
a variety of quantum systems such as in nuclear motion 
in diatomic molecules, in Josephson junction devices, and 
in active single-electron atoms. 

Let us concentrate on the latter example in the form 
of the simplest single electron atom, i.e., the hydrogen 
atom. Due to the application of a linearly polarized 
laser field F(t), the electron dynamics is effectively con- 
fined to one dimension along the polarization axis. The 
Hamiltonian for such a simplified (yet reliable de- 
scription of the hydrogen atom, which is here also per- 
turbed by a stochastic force £(i), reads as (atomic units, 
h = m = e = 1, are used unless stated otherwise) 
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H(x,t) = ^ + V(x) + x{F(t)+ffl}, (1) 

where x is the position of the electron and p = 
—i d/dx is the momentum operator. The potential 
is approximated by a non-singular Coulomb-like form 
V{x) = — \j\Jx 1 + a 2 . Such a soft-core potential with 
parameter a has been routinely employed to study atomic 
dynamics in strong laser fields [15|. It successfully 
describes many experimental features of multiphoton 
or tunnel ionization 14j, and the observation of the 
plateau in the higher harmonic generation spectra [l5j |. 
The external perturbations (last term in Eq. (1)) is 
dipole-coupled to the atom. The laser field is a non- 
resonant femtosecond pulse (duration 20 optical periods) 
described as, F(t) — f(i)Fo sm(ujt + 5). Here f(t) defines 
a smooth pulse envelope with Fq and lo denoting peak 
amplitude and angular frequency, respectively. The noise 
term £(t) is a zero-mean (£(t)) = 0, Gaussian white noise 
with autocorrelation function 
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and intensity D [lij ]. 
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Due to the stochastic nature of the Hamiltonian, the 
quantum evolution is nondeterministic. Thus an averag- 
ing over a large number of realizations of the stochas- 
tic force is required in order to produce a statistically 
meaningful solution of the time-dependent stochastic 
Schrodinger equation 
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H(x,t) #(x,t). 
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For a given realization, the numerical solution 
of the Schrodinger equation amounts to propa- 
gating the initial wave function \^q) using the 
infinitesimal short-time stochastic propagator, 

U^(At) — exp ^— i Jt +At H (x, t)dt\ . One can com- 
pute Up (At) using the standard FFT split-operator 
algorithm [It}, with the stochastic integration in the 
exponential interpreted in the Stratonovitch sense 
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Successive applications of the propagator U^(At) 
advance l^o) forward in time. Note that the initial 
state | Wo) is the ground state of the system having 
a binding energy of lb = —0.5 a.u.. This is obtained 
by the imaginary-time relaxation method for a 2 = 2 
14j . To avoid parasitic reflections of the wavefunc- 
tion from the grid boundary, we employ an absorbing 
boundary 17]. The observable, such as the ionization 
flux leaking in the continuum on one side, is defined as 
jR{xR,t) = i?e[W* p 9] XR , where xr is a distant point 
(typically 500 a.u.) near the absorbing boundary. The 
ionization rate is integrated over a sufficiently long time 
interval to obtain the corresponding total ionization 
probability, P = J R (x R ,t)dt. 

First, we will discuss the response of the atom interact- 
ing with a short but strong laser pulse only. It produces 
(nonlinear) ionization of the atom which is most easily 
understood, especially in the time domain, with the pic- 
ture of tunneling ionization. Ionization flux is produced 
close to those times when the potential is maximally bent 
down by the dipole-coupled laser field. This is illustrated 
in Fig. 1 with the temporal evolution of the ionization 
flux for a 20 cycle laser pulse (shown in the top parts of 
Fig. 1) with two different peak amplitudes F = 0.05 a.u. 
[Fig. 1(a)] and F = 0.02 a.u. [Fig. 1(b)]. Time-resolved 
ionization peaks separated by the optical period (2it/uj) 
are clearly visible for both peak field amplitudes. In ad- 
dition, Jn(t) shows a complex interference pattern due 
to the modulated Coulomb barrier in Fig. 1(a). How- 
ever, quite strikingly, if F is reduced to 0.02 a.u., the 
ionization collapses by around five orders of magnitude 
as shown in Fig. 1(b). One can conclude that the pho- 
toionization dynamics is highly nonlinear, and in partic- 
ular it exhibits a form of "threshold" dynamics where 
the threshold is created by the condition for over barrier 
ionization. 

Almost nothing is known about the quantum ana- 
log of SR in such "threshold-crossing" quantum devices 
[2|. Naturally, in the context of SR the question arises 
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FIG. 1: Nonlinear ionization flux Jn(t) (lower part of plots) 
induced by a 20 optical cycles laser pulse shown in the respec- 
tive upper plots with ui = 0.057, 8 = and peak amplitude, 
F = 0.05 (a), and F = 0.02 (b). The pulse envelope f(t) 
is unity except single-cycle sinusoidal rising and falling edges. 
The threshold for over barrier ionization is Fth = 0.067 a.u.. 



if noise can recover the strong periodic ionization flux 
for the "weak" laser pulse. To answer this question, we 
show the ionization flux averaged over 50 realizations in 
Fig. 2(a) when a small amount of noise is added to the 
weak laser pulse (F = 0.02). Note that the noise is 
switched on for the same time interval as the 20 cycle 
laser pulse. One can see that for very small noise am- 
plitude \pD — 0.00024, the periodic structure in atomic 
ionization gets amplified by one order of magnitude as 
compared to the case of the laser pulse alone [compare 
Fig. 1(b) with Fig. 2(a)]. However, one might ask if 
such an enhancement could be due to the noise alone. 
Contrary to the coherent excitation, the noise alone pro- 
duces an, on average, featureless ionization profile [see 
Fig. 2] . Here the noise causes ionization from the atomic 
ground state, which is different from the studied inization 
of the Rydberg atoms [l8(. This purely stochastic ion- 
ization for the feeble noise \[~D = 0.00024 is considerably 
smaller, than the corresponding case of laser pulse with 
noise [Fig. 2(a)]. Hence, the observed net enhancement 
can be attributed to a nonlinear quantum interaction be- 
tween coherent pulse and noise, which is also beyond the 
response to a simple quantum addition of the individual 
external fields. 

As the noise level is further increased, we observe an 
enhancement of the periodic ionization profile by many 
orders of magnitude as shown in Fig. 2(b)-(c). However, 
the increase in noise level also causes the stochastic ion- 
ization curve to rise rapidly. Eventually, for strong noise 
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FIG. 2: Ionisation flux for a weak laser pulse Fo = 0.02, with 
three values of noise amplitude \[D (a) 0.00024, (b) 0.0015, 
and (c) 0.018. Background featureless curves show the cor- 
responding purely noise-driven (Fo = 0) flux. The flux is 
averaged over 50 realizations. 



have plotted the enhancement factor rj versus the noise 
amplitude. It exhibits a sharp rise, followed by a maxi- 
mum at a certain value of the noise (point B) , and then a 
gradual fall off. It is worth mentioning that only a mod- 
est noise to laser ratio (^D opt /F — 0.075) is required 
to reach the optimum enhancement (here rj max — 36), as 
indicated by a typical optimal pulse shape in the inset of 
Fig. 3. 

The enhancement in photoionization can be under- 
stood by a simple two-step mechanism. First, starting 
from the ground state, the atom absorbs energy from 
the noise leading to a Boltzman-like population of energy 
levels. In a second step, the laser field causes ionization 
from the electron wavefunction Boltzmann distributed 
over many (excited) states. This simple picture is in- 
deed verified by separating the laser irradiation and the 
noise input in time. A sequential application of the noise 
followed by the laser pulse leaves the SR curve almost in- 
variant. On the contrary, by reversing the sequence, i.e., 
laser pulse first and then noise, destroys the SR. Note 
that due to the level structure of the atom (many bound 
states and continuum) there is no relevant noise induced 
time-scale as known from two-level systems. Thus, the 
optimum enhancement does not show the characteris- 
tic synchronization between coherent and noise induced 
time-scales as in the TLS [H, Q ■ Indeed, in some classical 
systems SR has been shown to exist without any explicit 
synchronization requirement [X9j | . Atomic ionization un- 
der a driving laser field provides the quantum analog to 
synchronization free SR. 

It is worth mentioning that the features presented here 
are robust with respect to the choice of parameters. For 
example, we have observed the quantum SR for pulses 
lasting from a few cycles to few hundred cycles, and for 



case, the coherent structures tend to wash out and the 
noise dominates the ionization [Fig. 2(c)]. This suggests 
the existence of an optimum ratio between the noise and 
laser amplitudes which leads to a maximum ionization 
enhancement. 

To quantify the quantum stochastic enhancement we 
define the enhancement factor 
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with Pq = P s +P n . Here P s + n denotes the average ioniza- 
tion probability (IP) due to the presence of the laser pulse 
with noise, P s and P n are the individual IP for laser pulse 
and noise, respectively. Although this is different com- 
pared to the quantifiers commonly used 0, 0| , rj is more 
suitable for our case. One can verify that a zero value 
of rj corresponds to the case when either the laser pulse 
(P s ^> P n ) or the noise (P s <C P n ) dominates. Further- 
more, rj characterizes a truly nonlinear quantum inter- 
play as it also vanishes if we assume a "linear" response 
as a sum of individual IP, P s + n = P s + P n . In Fig. 3 we 
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FIG. 3: The enhancement in photoionization due to quantum 
SR. The points marked A-C correspond to the noise ampli- 
tudes of Fig. 2(a)-(c), respectively. Inset: A typical optimal 
pulse for point B with ^/D^ = 0.0015, F = 0.02. 
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and molecular systems, our results might provide valu- 
able insight into the possible role of noise in designing op 



timal quasi-coherent quantum control schemes llj, [21 1 . 
Finally, analogous effects are also expected in other sys- 
tems such as in photo-fragmentation of anharmonic di- 
atomic molecules 22 , [23| , and in the recently observed 



multiphoton transitions in current-biased Josephson de- 
vices [241 1 . provided a coupling with an incoherent per- 
turbation exists. 

We thank A. Kenfack, W. Peijie, N. Singh, A. Buch- 
leitner, and P. Hanggi for fruitful discussions. 



FIG. 4: Enhancement induced by a broadband chaotic light. 
The peak amplitude and frequency of the 20 cycle laser pulse 
are Fo = 0.02 and lu = 0.057, respectively. The bandwidth 
of chaotic light is Au) — 0.75 with central frequency u>o = 
0.375. Inset: power spectral density (PSD) of the chaotic 
light compared to the one for the white noise. 



lo ranging over more than one order of magnitude from 
infrared to near UV frequencies. 

The intensity required for the white noise renders 
an experimental test of our prediction quite challeng- 
ing. However, one could think of replacing the white 
noise by a quasi- white noise, i.e., chaotic light with a fi- 
nite but broad bandwidth Aw. This can be realized by 
adding a large number (here 1024) of independent phase- 
randomized frequency modes. The inset of Fig. 4 shows 
an example of a characteristic spectral density of such a 
broadband source having Aw = 20 eV (0.75 a.u.). Mod- 
ern pulse shaping techniques can generate chaotic light 
pulses with a bandwidth as large as 30 eV [2(| Hl[ . In 
Fig. 4 we have plotted r\ defined in Eq. (4) vs the rms 
amplitude of the chaotic light (analogous to the noise 
amplitude) for the previously used coherent pulse. One 
can see that such chaotic light can preserve the features 
of quantum SR, with an almost identical optimum com- 
pared to the one for the white noise case of Fig. 3. This 
observation may creat new possibilities in quasi-coherent 
control schemes of similar quantum systems. 

In conclusion, we have demonstrated a new form of 
quantum stochastic resonance in the dynamics of the 
simplest atomic system for the first time. This gener- 
alized quantum SR leads to a dramatic enhancement (by 
several orders of magnitude) in the nonlinear ionization 
when a modest amount of optimum white noise is added 
to the weak few cycle laser pulse. The same effect is also 
achieved if one uses (realizable) broadband chaotic light 
instead of white noise. We emphasize that the effet is 
robust with respect to a range of experimentally accessi- 
ble parameters. In addition to substantially broadening 
the existing paradigm for quantum SR to generic atomic 
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